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C^ I Abstract. We describe a simple derivation of virial relations, for arbitrary physical systems 

that are governed by an action. The virial theorem may be derived directly from the action. A, 
with no need to go via the equations of motion, and is simply a statement of the stationarity 
of the action with respect to certain variations in the degrees of freedom. Only a sub-class 
of solutions that obey appropriate boundary conditions satisfy each virial relation. There 

is a set of basic virial relations (one for each degree of freedom) of the form -g^|^_^ — 0, 

/^ ' where Al is the action in which one degree of freedom has been scaled by a constant factor c. 

^— N ' Linear combinations of these may be put in simple forms, taking advantage of homogeneity 

- ' properties of the action, and of dimensional considerations. When some of the degrees of 

freedom are of the same type a tensor virial theorem presents itself; it may be obtained in 

a similar way by considering the variation of the action under linear transformations among 

these degrees of freedom. Further generalizations are discussed. Symmetries of the action may 

lead to identities involving the virial relations. Beside pointing to a unified provenance of the 

virial relations, and affording general systematics of them, our method is a simple prescription 

,S^ • for deriving such relations. It is particularly useful for treating high-derivative and non-local 

CIh' theories. We bring several examples to show that indeed the usual virial relations are obtained 

' I by this procedure, and also to produce some new virial relations. 

^ ' 1. Introduction 

C^ ■ 

In default of a general definition of the virial theorem we describe it, drawing from known 
examples (see, for some of many, [1-9], and, in particular, [10]): It consists of a set of global 
k><( , (integral) relations that are satisfied by a subclass of solutions of the equations of motion. 

'Oj . This sub-class is defined by requiring that the solutions obey certain boundary requirements, 

j^ ' The derivation of these relations proceeds as follows: One contracts the equations of motion 

with functions of the degrees of freedom and integrate over the variables on which the latter 
depend. One then integrates by parts so as to reduce the order of derivatives, as many times 
as possible-discarding boundary terms, as one proceeds, by imposing requirements on the 
boundary behavior of the solutions. One ends up with a virial theorem consisting of (1) the 
set of integral relations, and (2) the set of boundary conditions under which they apply. 

The applications of this poor-man's substitute for the equations of motion are many and 
varied (see e.g. [3][10][11-13], and the many uses described in [10] ). Its usefulness draws 
partly from the fact that it involves derivatives of the degrees of freedom of lower order than 
appear in the equations of motion. (We shall show, in fact, that the order of derivation 
appearing is the same as that in the action itself.) Thus, for example, the applications in 
astrophysics require knowledge of only particle velocities, which can be measured, and not of 
the accelerations, which cannot. This feature also makes the virial theorem a useful tool for 
checking numerical solutions of the equations of motion. 

Here, we give a systematic treatment of the virial theorem, for arbitrary systems governed 
by equations that are derived from a variational principle, in a way that highlights its origin, 
and greatly facilitate its derivation in the case of higher order theories, and especially for 
non-local theories. 

We explain how the virial relations emerge from the action in sections 2 and 3, showing 
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that they are tantamount to stationarity of the action under certain increments in the degrees 
of freedom; when the variation of the action under such increments may be expressed simply, 
a convenient virial relation results. We demonstrate the procedure with various examples in 
section 4. 



2. Derivation of the virial relations from the action 

Virial relations may be derived for action-governed systems as follows: The system is 
described by d degrees of freedom /i['7*''], 1 < ^ < c?, where, for each i, rj'-^^ is the set of 
variables , rji , ..., on which /; depends (time, path length, space-time coordinates for fields, 
Fourier variables if the degree of freedom is described in some Fourier space, etc.). One 
assumes that there exits an action A: a functional of fi , such that under a variation Sfi in fi 



5A = Y,jdr]^'^ Q^ ^/, + ^[^,('^/,)]e„d., 



(1) 



where Qi{fi, ..., fdiV^''^) is a functional of the degrees of freedom and a function of 77''' only; 
J^i is a linear operator acting on 5fi, and may depend as a functional on all the /^s. The 
term [ ]ends is a boundary term evaluated at the boundary of the respective 77^*' space. One 
then postulates that the physical states are described by /^s for which A is stationary under 
variations 5fi that annihilate the ends term. These /^s are then the solutions of the equations 
of motion 

g,(/i,...,/d,77«) = l<i<d. (2) 

Consider now, as a specific change, a rescaling of a single degree of freedom /j 

h^{l + e),U = ,h + 5,U. (3) 

The change in the action is then, by eq. (1) 



5A^e / ci7?W Q,h + e[T,{h)Uds- 



On the other hand, we can write 



5A^e 



dA: 



where 

We thus have 



dA 



dc 



) ^Ji: • ■• 1 Id) • 



dc 



drj 



{i) 



.h + [T^{h)Uds. 



(4) 

(5) 
(6) 
(7) 



Equation (7) is but an identity. If we now specialize to solutions of the equations of 
motion (2) we get 

dAi 

— [^i{fi)]ends, (8) 

c=l 



dc 



which is non-trivial (i.e. is not satisfied by generic /J??*''] that are not solutions of the 
equations of motion ). It tells us that certain expressions derived from the action are pure 
boundary terms. 

The virial relation results when we further restrict ourselves to solutions for which the 
ends terms may be discarded. This may happen because the particular /^s satisiy boundary 
conditions in 77^'' space that render [Ti{fi)]e„ds zero, or, more generally, because for a large 
enough ry^'' volume the ends term remains finite while A becomes infinite (in which case we 
normalize by the volume, and the ends term vanish in the limit). We then obtain as our 
desired virial relations 

BA^ 
^ =0 l<i<d, (9) 

c=l 



dc 



2 



each satisfied on its own sub-class of solutions for which the ends term can be discarded. 
Clearly, cq. (9) is just the relation we would get by taking the equation of motion derived by 
varying A with respect to /;, then multiplying by fi, then integrating by parts as many times 
as needed to reduce the order of derivatives back to that in A. The present method is then 
but a way to shortcut this return journey to the equations of motion. 

Linear combination of relations (9) with coefficients a^-which should hold under the 
appropriate boundary behavior-may be conveniently written as 



oa: 



do 



- 0, (10) 



c=l 



where 

^:i-"<'=^(c"i.A,...,c"d/rf). (11) 

Some of these linear combinations may be more useful than others. For example, some may 
be reduced to simpler forms by employing homogeneity properties of the action, or its parts. 
If there is a choice of a^ for which A"'^"'°''^ = c^A, with /3 7^ 0, then, by eq. (10) , ^ = is 
one of the virial relations. (If /3 = the d relations (9) are dependent, but the theory is then 
degenerate.) This is the case for any system with a Newtonian kinetic action and harmonic 
forces, in which case ^ = implies the equality of mean kinetic and potential energies. It is 
also the case for the Einstein-Hilbert action in vacuum, in which case the resulting relation 
is just the integral of the trace of the equation of motion. 
More generally, the action may be split into terms 

a 

all homogeneous with the same set of ct^s, but each with its own (3a. Then 

^/3.A = (13) 

a 

is one of the virial relations. This is the case for a system of particles with a Newtonian kinetic 
action, and forces that are homogeneous in the positions (e.g power-law forces). Equation 
(13) then tell us that the ratio of mean kinetic and potential energies is fixed for the theory. 
It is also the case for General Relativity with matter in the form of charged particles, and 
electromagnetic fields (see sec. 6). 

We may also apply dimensional analysis to obtain some useful combinations of the virial 
relations. Suppose the action depends on some constants g^, 1 < p < n. Consider then, 
for instance, how the different quantities change under a change in the units of length by a 
factor c. Depending on their respective dimensions, fi will change to c"^ fi, Qp will change to 
c^'fjfp, while A will change to c'^A. (We assume that 77^'' have no length dimension, as when 
they denote the time, or if they do we redefine them not to have such dimension, and then 
we have to add to the qp another constant of the dimensions of length.) Dimensional analysis 
then tells us that 

^(c"i/i,-,c^^'Zi,-)=cM, (14) 

and so 

^.i...«, ^ ^f3^(^f^^ ^ j^^ ^--^^^^ ^ c-^"(?„). (15) 

From this and eq. (10) we then obtain 

dA 
/3^-^7p9p7^=0, (16) 

OQp 

as one of the virial relations. This is particularly useful when A depends only on a small 
number of constants. 



3. Tensor virial theorems and further generaUzations 
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Often, some of the degrees of freedom are of the same type-say fi for 1 < i < K, 
collectively designated /. Then, further virial relations suggest themselves. (Examples of 
degrees of freedom of the same type are the position coordinates of a particle, or indeed, of 
different particles, the components of the electromagnetic vector potential field, the elements 
of the metric tensor etc..) 

We now consider the variation of the action under variations of the form 

f^ Mf, (17) 

where M = I + £ (with £ infinitesimal) is a constant K x K matrix. From eq. (1) we have 



SA = Y, ^^j J dv 0. /. + [^^if,)l 



(18) 



(All the Qi and fi, i < K, arc defined on the same rj space.) 
On the other hand, defining 

^m(/,/k+i,...,/d) ^^(M./,/k+i,...,./d), (19) 

we see that 



M = ^£.,y,,, (20) 



where 



V,. 



dA 



M 



dM.j 



(21) 



Thus 

V,, = JdrjQ, /,. + [Td.f,)U,s- (22) 

As in section 2, we confine ourselves to solutions of the equations of motion for which the 
ends term in eq. (22) can be discarded and then 

V,, - 0, (23) 

which is the desired K x ii'-tensor virial thcorcm(thc diagonal elements of relation (23) are 
the same as eq. (9) ). In addition wc still have, of course, the relations of type (9) for all the 
degrees of freedom fi i > K. 

Symmetries of the action under transformations of / lead to constraints on the tensor 
virial relations. For example, say that A is invariant under a rotations in / space, i.e., under 
/ — > Uf, for any orthogonal matrix (continuously connected to the unit matrix). Then, from 
eq. (20) , J2i-^ij^ij = ^'^^ 9,ny antisymmetric matrix £ (the generators of the orthogonal 
transformations). Thus, V must be a symmetric matrix: Vij = Vji is then an identity, and 
there are only K{K + l)/2 independent relations. Note that symmetries of the action do not 
necessarily imply symmetries of the above type where only some of the degrees of freedom 
are transformed. 

The information held by the equations of motion on their solutions is infinitely greater 
than that supplied by the finite number of discrete, global relations that we have discussed so 
far. It should be obvious then that one can derive an infinite number of further such relations. 
These emerge if we consider the variation of the action under general infinitesimal increments 
of the form 

/. ^/.+6(./i,-,/d). (24) 

Here ^i may be a function of the degrees of freedom and also a functional of them, but it must, 
of course, be "of the same type" as fi itself (in the sense that the action is defined if we replace 
fi by ^i). As before, if the ends term [J-'i{S,i)]ends rnay be discarded, the variation of A under 
(24) vanishes on solutions of the equations of motion. When, for given A and £,, SA can be 
written in a simple form-as was the case for the linear increments we had discussed earlier-a 
useful relation may be obtained. The point is that because the increment (24) is defined in 
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terms of the degrees of freedom themselves it enables us to formulate the extrcmum condition 
on A, and the required boundary conditions, in terms of the /^s themselves thus leading to 
a virial-like relation. When it is not necessary to have a relations written in terms of system 
properties alone, it may be useful to consider more general increments in which ^^ is also an 
explicit function of t]''-*. Thus considering fi -^ ^ijkfjfk will lead to higher-tensorial-order 
relations (for a system of Newtonian particles they correspond to relations connected with 
higher moments of the mass distribution-see section 5). 

While these further relations are not as amenable to many of the uses as their predeces- 
sors, they may be useful, depending on the particular case. For example, in connection with 
checks on numerical solutions of the equations of motion the may provide checks of aspects 
not covered by the basic relations (see section 4). 

We have discussed virial relations as they apply to the full sub-class of solutions that 
satisfy the appropriate boundary conditions. When one limits himself further to special 
solution-such as ones with certain symmetry properties-the virial theorem may, in general, 
be simplified; we shall not deal with such further restrictions. 

4. Examples 

Now let us look at some examples, mainly to demonstrate the procedure and various 
relevant points discussed generally in the previous sections. We start with the archetype 
of all virial theorems: that for a system of N Newtonian particles interacting via a general 
iV-body potential. The action for a finite time span (ti, ^2) is 

A = Ak + Ap, (25) 

with the kinetic action 

A.=T-^iy^Y.-^.[f)\t, (26) 

and the potential action 

Ap = -r-i/'$[ri(f),...,f^(t)] dt. (27) 

We have normalized the standard action hy T = t2 — ti, to make A finite for T -^ 00. 

The degrees of freedom are the 3N components rf of the positions. We get, straightfor- 
wardly, from eq. (9) the 3N virial relations 

(m,(ff )2) - (rfa$/9rf) = 0, (28) 

where 

(Q) =T-' r Q dt. (29) 



Usually this average is taken in the limit t-i -^ — cxd, and ^2 -^ 00, and it then holds for 
solutions in which all particles are bounded in a finite volume. The standard scalar virial 
theorem for this case: 

(^m,«,2)_(^f^.V,^$)=0, (30) 

i i 

is obtained by differentiating ^[cri(i), ..., cr]v(i)] with respect to c (at c = 1). There is a 
?>N X 3iV tensor virial theorem 

{m,fy^^) - (rf 9$/arf) = 0, (31) 

which is obtained from eq. (23) , by considering r" ^r r^ -\- er^, for fixed values of all indices. 
Neither term is symmetric when the masses are not all equal, because neither term in the 
action is symmetric under rotations in configuration space. The standard tensor virial theorem 



is obtained by taking the trace over the particle indices, or by considering r" ^ r° + er^ for 
aU values of i but fixed a, (3. 

{ E "^'^r^O ^ i^td^^/drf) - 0. (32) 

i 

The kinetic term is now identically symmetric because the kinetic action is symmetric under 
simultaneous rotations of all the Vi. The antisymmetric part of the potential term is ^^ r^ x 

Vp.^ which vanishes identically if $ is rotationally invariant, in which case relation (32) is 
identically symmetric in a/3. 

We now add a field degree of freedom by considering a system of N gravitation particles 
of masses m; and we also want to solve for the (Newtonian) gravitational field (p{r) which is 
thus also a degree of freedom. The action is 



I, i i 



[Vififd'r \ . (33) 



It gives the standard equation of motion of the particles in the potential (p, and variation over 
Lp gives the Poisson equation for ip with the density p{r) = ^^ mi5^{f — Ti). 

Beside the scalar and tensor virial relations we wrote in the previous example (with 
$[ri(t), ...,rjv(0] = Si "^i'/'i^ilO]); '^6 no"^ have a relation resulting from p ^ cp which is 
read directly from the action-the three terms in A being homogeneous in p, with powers 0,2, 
and 1 respectively: 

{Y^mMn)) + ^{j{ypfd^r)^Q. (34) 

i 

To demonstrate the potential usefulness of relations resulting from the more general types of 
variations eq. (24) , consider an infinitesimal increment 

p ^ eJ{p). (35) 

Substituting in the action, and taking the first order in e we get 

( E m,J[p{n)]) + ^{j J'{p){ypf<i'r) = 0, (36) 

i 

with the boundary requirement being that J{p) = 0. Clearly, this infinite set of relations 
captures more of the information in the equation of motion. By choosing J properly we may 
accentuate different regions of space, for example regions where there are better measurements, 
or where better approximations can be made. 

Consider now a system of particles governed by a general kinetic action, so the action is 
of the form 

A^AK[ri{t),...,fN{t),qi,...,q„]-T-' f <i>[r,{t), ...,rr,{t)] dt. (37) 

where Ak is a general functional of the particle trajectories ri{t) (including higher-derivative 
or, indeed, non-local actions, in which case the full world lines of the particles are integrated 
on), and Qp are constants. We apply eq. (16) noting that j3 = 2 (with our choise of normal- 
ization A has dimensions of energy), and that $ must be of the form 

$-$„(/.(fi//,...,nv//). (38) 

Here, $o, and / are constants with the dimensions of energy and length, respectively. Thus, 
we obtain in generalization of the scalar virial theorem resulting from ri{t) -^ (1 + e)ri{t) (see 
[13]) 

dA ^ 

2^^ - E 7p?p Trp - ( E ""^ • v., $) = 0. (39) 



As in eq. (16) , 7p are the powers of length in the dimensions of the constants q^. 
When Ak is a Lagrangian action: 

Ak=T-^ I Lj,[r,y^\7f\...,q;\ dt, (40) 

with the kinetic Lagrangian L^ depending on higher time derivatives 77 of the trajectories, 
we can write, in generahzation of the standard scalar virial theorem, either 

(EE^t'-i^)-(E^>v,a>).o, (41) 

obtained from the variations of A under r* (t) -^ cfi{t), 1 < i < N, or, equivalently, from eq. 
(39), 

" Pi/ T \ ^ 

2{L^) - J2 7p9p-^ - ( E ^^^ • ^^^'^) = 0- (42) 

The latter is particularly useful when L^ depends only on a small number of constants. 
The tensor virial theorem, also read directly from A is 

^--(EE^r«|^)-(E^"S) = o (43) 

i— 1 I ^ i — 1 ^ 

(r° is the ^th time derivative of the a component of r^). Rotational invariance of Ak implies 
that the kinetic part of Vafj is symmetric, as an identity, although this is not manifest in eq. 

(43) • 

For a general relativistic system of particles of masses rrii and electric charges e^, the 

action is [14] 

A = Ag + Ap + Aem + A,, (44) 

where 

is the Einstein action for the metric; 

1/2 



dx^dxl 
"■'''"' ^""'^ dX d\ 



(46) 



Ap = — y^ nii / dX 

i 

is the particle action; 

Aem = -\ J {A,,, - A^,.)3''"5""'(^;3.c ~ A„,^)g^/'d*x (47) 

is the action for the electromagnetic field; and 

A-^e.y dA^A^(x,) (48) 

i 

is the field-particle-interaction contribution. The degrees of freedom here are the particle 
space-time coordinates a;f (A), the four electromagnetic-potential fields A^(x), and the com- 
ponents, g^i,(x), of the metric tensor. 

Under g^^, -^ cg^^ the terms in the action are homogeneous: Ac -^ cAq, Ap -^ c^^'^Ap, 
and Aem, and Ae are invariant. Thus we get as one of the virial relations 

Ag + \Ap = 0. (49) 



This can be shown to be just the space-time integral of the trace of the (Einstein) equation 
of motion. This result is easily seen to be general: The definition of the energy-momentum 
tensor, T^^, is such that under a change Sg^^ we have 5A = (1/2) J Sg^^j^T^'^g^^^d'^x [14] ; so, 
under g^^, ^ (1 + e)g^^, SA = {e/2) J Tg^^^d'^x, for whatever the source is. The boundary 
requirements for this type of relation is satisfied if the metric becomes flat at spatial infinity 
(so no gravitational waves escape to infinity), and if there are periodic boundary conditions 
in time (e.g. a stationary solution). 

There is also homogeneity under A^, -^ cA^, which gives a virial relation 

2Aem + a = (50) 

(requiring that no EM waves escape to infinity). 

There are tensor virial theorems for (7^^ and for A^,; the latter, for example, gotten by 
considering Aa -^ Aa + eAp reads 

/r fix" 
A/3,m5'*"5"''(A.. - A,,^)g'/'d'x + E ^^ 7 dX^A,{x,) - 0. (51) 

i 

All these relations may be easily derived from the equations of motion, as is always the case 
with low-order theories. We bring them here to show how they follow from our systematic 
procedure, by mere inspection of the action, and without having to derive the equations of 
motion first. 

By employing a more general conformal increment 

eq. (49) generalizes to 

f[-{STrG)-''R + T]Jg^^''d*x^O, (53) 

and there is a similar generalization for eq. (50) . 

The action may have terms quadratic in the curvature tensor-as is the case for the first 
order semi-classical corrections to General Relativity obtained from string theory-such as 

f R^g'^^d^x, f R.^R^-'g'^^d^x, j R^^.^R'^^^'^g^'^d^x. (54) 

These are invariant under g^^ — > cg^^, and do not modify relation (49) . In order to capture 
them in a virial relation we may need to employ a tensor relation, or use one that is induced 
by increments of type (24) . 
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